Charmonium mass in nuclear matter 
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The mass shift of charmonium states in nuclear matter is studied in the perturbative QCD 
approach. The leading-order effect due to the change of gluon condensate in nuclear matter is 
evaluated using the leading-order QCD formula, while the higher-twist effect due to the partial 
restoration of chiral symmetry is estimated using a hadronic model. We find that while the mass 
of J/ip in nuclear matter decreases only slightly, those of ^'(3686) and i/)(3770) states are reduced 
appreciably. Experimental study of the mass shift of charmonium states in nuclear matter can thus 
provide valuable information on the changes of the QCD vacuum in nuclear medium. 
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Understanding hadron mass changes in nuclear 
medium and/or at finite temperature can provide valu- 
able information about the QCD vacuum [1-3]. It is also 
relevant phenomenologically to the interpretation of ex- 
perimental results from relativistic heavy ion collisions 
[4], in which a hot dense matter is formed during the 
collisions. Previous studies have been largely concerned 
with hadrons made of light quarks [3]. Only recently 
were there studies of the in-medium masses of hadrons 
consisting of heavy charm quarks. Using either QCD sum 
rules [5,6] or the quark-meson coupling model [7], it has 
been found that the mass of D meson, which is made 
of a charm quark and a light quark, is reduced signifi- 
cantly in nuclear medium as a result of the decrease of 
the light quark condensate. For the J/V', which consists 
of a charm and anticharm quark pair, both the QCD sum 
rules analysis [8] and the LO perturbative QCD calcula- 
tion [9,10] show that its mass is reduced slightly in the 
nuclear matter mainly due to the reduction of the gluon 
condensate in nuclear matter. 

The change of hadron masses at finite temperature is 
best studied using the lattice gauge theory as its predic- 
tion is less model dependent. Recent lattice gauge cal- 
culations at finite temperature with dynamical quarks 
have shown that even below critical temperature the 
interquark potential at large separation approaches an 
asymptotic value Voo{T) that decreases with increasing 
temperature [11]. This transition from a linearly rising 
interquark potential in free space to a saturated one at 
finite temperature is due to the decrease in the string 
tension and the formation of Qq and qQ pairs, where q 
denotes a light quark, when the separation of the two 
heavy quarks (Q) becomes large. The decrease in Voo{T) 
can thus be interpreted as a decrease of the open heavy 
quark meson (Qq or qQ) mass mjj, such as the D meson 
mass, at finite temperature [12,13]. Furthermore, the de- 
crease of rriH seems to be a consequence of the reduction 
in the constituent mass of light quark as the tempera- 
ture dependence of Voo{T) is similar to that of the chi- 
ral condensate (qq) [14]. This relation between the mass 



rriH and the chiral order parameter also follows naturally 
from the heavy quark symmetry [15]. With the finite 
temperature interquark potential, it has been shown via 
the solution of corresponding Schrodinger equation that 
the masses of charmonium states are reduced at finite 
temperature as well [16,17]. 

At finite density, lattice gauge calculations are at 
present not feasible for studying the heavy quark poten- 
tial or the mass mn of the open heavy quark meson. 
Masses of these heavy quark systems are, however, ex- 
pected to change appreciably in nuclear medium. Model 
independent estimates have shown [18,19] that conden- 
sates of the lowest dimensional operators (^G^) and (qq) 
decrease, respectively, by 6% and 30% at normal nuclear 
matter, which are significant changes expected only near 
the critical point at finite temperature [20] . As in the case 
of finite temperature, the reduction of gluon condensate 
leads to the softening of the confining part of interquark 
potential [21], while the decrease of quark condensates 
implies a drop of the open heavy quark meson mass or 
Voo of the heavy quark potential. Both are expected to 
lead to nontrivial changes in the binding energies of the 
charmonium states ■0(3686) and V'(3770), as their wave 
functions are sensitive to both the confining part and the 
asymptotic value of the interquark potential. 

In this letter, we evaluate the mass shift of ?A(3686) 
and ■0(3770) due to changes in the gluon and quark con- 
densates in nuclear medium. The effect of the gluon con- 
densate is determined using the leading-order LO QCD 
formula, which was developed in Refs. [9,22] and has been 
used to study the J/^p mass in medium [10]. The effect 
due to the change in quark condensates is difficult to cal- 
culate using the quark and gluon degrees of freedom as 
they appear as higher twist effects in the operator prod- 
uct expansion [23,24]. However, its dominant effect to a 
heavy quark system is to reduce Voo as a result of the 
decrease of D meson in-medium mass, about 50 MeV in 
normal nuclear matter due to the 30% reduction in the 
quark condensate [5-7,25]. Therefore, we can study the 
effect of changing quark condensate on the charmonium 
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states at finite density by using a hadronic model to cal- 
culate their mass shifts due to the change of D meson 
in-medium mass. Combining the effects from changes 
in the gluon condensate and in mo, wc find that both 
^(3686) and -0(3770) masses are reduced appreciably at 
normal nuclear matter density. 

The mass shift of charmonium states in nuclear 
medium can be evaluated in the perturbative QCD when 
the charm quark mass is large, i.e., rric — > oo. In this 
limit, one can perform a systematic operator product ex- 
pansion (OPE) of the charm quark-antiquark current- 
current correlation function between the heavy bound 
states by taking the separation scale (/i) to be the binding 
energy of the charmonium [9,22,26]. The forward scat- 
tering matrix element of the charm quark bound state 
with a nucleon then has the following form: 
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Here, C„ is the Wilson coefficient evaluated with the 
charm quark bound state wave function and {On) jv is the 
nucleon expectation value of local operators of dimension 
n. 

For heavy quark systems, there are only two indepen- 
dent lowest dimension operators; the gluon condensate 
((^G^)) and the condensate of twist-2 ghion operator 
multiplied by {{^Gaij,G")). These operators can 
be rewritten in terms of the color electric and magnetic 
fields: {^E^) and {^B^). Since the Wilson coefficient 
for {^B"^) vanishes in the non-relativistic limit, the only 
contribution is thus proportional to {^E'^), similar to 
the usual second-order Stark effect. We shall thus calcu- 
late the mass shift of charmonium states due to change 
of the gluon condensate in nuclear medium by the QCD 
second-order Stark effect [10]. 

The mass shift of charmonium states to leading order 
in density is obtained by multiplying the leading term in 
Eq.(l), by the nuclear density pjv- This gives, 



Am^(e) 



■II 



dipik) 



dk 



k 



k'^/ruc + e 



-E^ 



Pn 



N 



2m 



N 



(2) 



In the above, m^^ and are the nucleon mass and the 
nuclear density, respectively; {^E^)n ^ 0.5 GeV^ is the 
nucleon expectation value of the color electric field and 
e = 2mc — m^. In Ref. [9], the LO mass shift formula 
was derived in the large charm quark mass limit. As 
a result, the wave function 'ijj(k) is Coulombic and the 
mass shift is expressed in terms of the Bohr radius uq 
and the binding energy eo = 2mc — mj/^. This might be 
a good approximation for J /ip but is not realistic for the 
excited charmonium states as Eq.(2) involves the deriva- 
tive of the wave function, which measures the dipole size 



of the system. We have thus rewritten in the above the 
LO formula for charmonium mass shift in terms of the 
QCD parameters = 0.84 and = 1.95, which are 
fixed by the energy splitting between J/iJj and '0(3686) 
in free space [9]. Furthermore, we take wave functions 
of the charmonium state to be Gaussian with the oscilla- 
tor constant /3 determined by their squared radii (r^) = 
0.472, 0.96^, and 1 fm^ for J/ip, ^'(3686), and V(3770), 
respectively, as obtained from the potential models [27]. 
This gives f3 = 0.52, 0.39, and 0.37 GeV if we assume 
that these charmonium states are in the 15, 2S, and ID 
states, respectively. Using these parameters, we find that 
the mass shifts at normal nuclear matter density obtained 
from the LO QCD formula Eq.(2) are -8, -100, and -140 
MeV for J/ip, ■0(3686), and V(3770), respectively. 

Although the higher twist effects on the charmonium 
masses are expected to be nontrivial, the result for J/ip is 
consistent with those from other non-perturbative QCD 
studies, such as the QCD sum rules [8,24] and the ef- 
fective potential model [28,29], which are all based on 
the dipole interactions between quarks in the charmo- 
nium and those in the nuclear matter. To go beyond the 
leading order in the OPE, we need to calculate the con- 
tributions from higher dimensional operators in Eq.(l), 
which include light quark operators. Explicit calculations 
from QCD sum rules for J/^ up to dimension 6 opera- 
tors [24] show that the effect due to change in the con- 
densates of light quark operators at dimension 6, which 
include {qTqqTq) and (qDGq) [23,24], is unimportant for 
the mass shift of J/tp- However, such calculation cannot 
be easily generalized to the excited charmonium states 
■0(3686) and ^0(3770), where the sum rules do not ex- 
ist even in the vacuum. On the other hand, the higher 
twist effects due to the light quark operators can be es- 
timated by considering the coupling of the charmonium 
to the DD states as in the potential model for charmo- 
nium states [27]. Therefore, instead of summing up the 
non-convergent contributions from the change in the light 
quark condensates in the OPE of Eq.(l), we estimate its 
contribution by evaluating the charmed meson one-loop 
effect on the mass of a charmonium with in-medium D 
meson mass predicted from the QCD sum rules [5,6] or 
the quark-meson coupling model [7] . 

Following the studies in Ref. [30] on p — tt interactions 
and in Ref. [31] on 6 — K interactions, we use the fol- 
lowing Lagrangian for interacting charmonium ^0 and D 
meson: 

C=\ {\D,T>\^ - ml\D\^) - \f,^F^^ + \ml^,r, (3) 

where F^^ = d^ijji, — Q^V/j the charmonium field 
strength, D^ = d^- i2g^DDi^ti, and D = (£)°,D+). 

The coupling constant g^oD can be determined using 
the 3P0 model [32] . In this model, the coupling constant 
is proportional to the overlap integral between the rel- 
ative quark wave functions of charmonium and the two 



2 



outgoing charmed mesons as well as to a coupling param- 
eter 7, which characterizes the probabihty of producing 
a Ught quark-antiquark pair in the ^Pq state. The result 
can be read off from Refs. [33,34] and is given by 



the D meson loop docs not contribute to the real part 
of the charmonium self energy in free space. The mass 
shift of the charmonium at finite density is then given by 
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where q is the three-momentum of D mesons in the ip 
rest frame and r = P/Pd with (3d (/?) being the os- 
cillator constant for D meson (tp) wave function. For 
the /3's, the same values will be used as in the LO 
QCD calculation. The values for 7 and 13d are taken 
to be 0.35 and 0.31 GeV, respectively, to reproduce 
both the decay width of ?/'(3770) to DD and the par- 
tial decay width of -0(4040) to DD, DD*. and D*D* 
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V'(3686), and 3(^^(1 - 5(i+2r^) /3g. 

Because of its momentum dependence, g^ooiq) takes 
into account the form factor at the ipDD vertex and al- 
lows also the coupling of the charmonium to off-shell D 
mesons. For i/)(3770), it can decay to DD in free space, 
and its on-shell coupling constant is 5^ (3770) dd (9 = 
(m^/4 - m|,)^/^) — 15.4. The coupling constants at 
g = are 15.3, 18.7, and 16.8 for J/tp, ^(3686), and 
ip{3770), respectively. The value for Jf-ip coupling to 
D mesons is slightly larger than that estimated by the 
vector meson dominance model [36,37] and by the QCD 
sum rules [38]. As the D meson momentum increases, its 
coupling constant to J/ip has a simple exponential fall 
off due to the 15 quark wave function of the J/tp. In 
contrast, the coupling constants of excited charmonium 
states V'(3686) and ^(3770) to D mesons fall off exponen- 
tially with the D meson momentum but vanish at certain 
q"^ as a result of the nodes in the 2S or ID wave functions 
of the excited charmonium states [33,35]. 

Similar to the method introduced in Ref. [31], we have 
used the above Lagrangian to evaluate the self-energy of 
a charmonium due to the D meson loop. After perform- 
ing the energy integral in the rest frame of the tp, i.e., 
k — (m^,0), the invariant part of the polarization n(fc) 
then has the following form: 
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where is the in-medium D meson mass and V de- 
notes that only the principle value of the integral is eval- 
uated. The subtracted term in the above equation is a 
renormalization constant which is determined by requir- 
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FIG. 1. Mass shifts of charmonium states J/^' (solid 
curve), ^(3686) (long dashed curve), and ■i/'(3770) (short 
dashed curve) as functions of D meson in- medium mass m*o. 
The circle indicates the expected mass shifts at normal nu- 
clear matter density. 

In Fig. 1, we show the mass shifts of charmonium 
states as functions of m|). It is seen that the mass shift 
of ip{'in{)) is negative for small negative shift of D meson 
mass but becomes positive when the D meson mass drop 
is large. In contrast, the mass shift of ^'(3686) is neg- 
ative for all negative mass shifts of the D meson. This 
difference can be understood from Eq.(5), where the in- 
tegral is a convolution of the form factor g^^D^iq^) with 
the terms in the square bracket, which are singular when 
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The integrand 



thus changes signs whenever the D meson momentum 
q passes through these singularities and finally becomes 
negative when q^ is larger than any of the singularities, 
which correspond to the energies of the virtual interme- 
diate D meson states. As in second-order perturbation 
theory, the contribution is attractive when the energy 
of the intermediate state is larger than the charmonium 
mass. However, the form factor decreases exponentially 
with q^ and can even be zero, the large negative contri- 
bution expected for a constant form factor is suppressed, 
leading thus to an increase of the -0(3770) mass when 
m_D — m|) > 10 MeV. On the other hand, the singu- 
larity of the integrand in Eq.(5) in the case of V'(3683) 
occurs only when 2m^ falls below its mass and there- 
fore has only a small positive contribution when q^ is 
very small, leading to a reduction of its mass for any D 
meson mass shift. For J/V', we find that its mass only 
increases slightly with dropping D meson mass and de- 
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when m*j-, = mo- This ensures that pends weakly on m^. For m_D 
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is the expected mass shift of D meson at normal nuclear 
matter density, the mass shift of J ftp is about 3 MeV. 
This result is consistent with that from the QCD sum 
rules [24] and is also expected from the potential model 
[27], where the J/ip wave function has only a small DD 
component. We note that the density dependence of the 
mass shifts of charmonium states, particularly the excited 
ones, is nonlinear if we use a linearly density-dependent 
D in- medium meson = niD — 50 p/po MeV in the de- 
nominator of Eq.(5). On the other hand, the mass shift 
obtained from the LO QCD formula in Eq.(2) depends 
linearly on the nuclear density. 

Adding the mass shift from the D meson loop effect 
to the result from the LO QCD calculation, we find that 
masses of charmonium states are changed by the follow- 
ing amount at normal nuclear matter density: 



Am j/^ = -8 -I- 3 MeV, 



Ato^(3686) = -100- 
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30 MeV, 
15 MeV, 



where the first number represents the shift from the LO 
QCD while the second number is from the D meson 
loop. The above results thus show that masses of excited 
charmonium states are reduced significantly in nuclear 
matter, largely due to the non-trivial decrease of the in- 
medium gluon condensate in the LO QCD formula for 
their masses. 

The mass shifts of both -0(3686) and "0(3770) in nuclear 
medium are large enough to be observed in experiments 
involving p — A annihilation as proposed in the future 
accelerator facility at the German Heavy Ion Accelerator 
Center (GSI) [39]. In these experiments, 0(3770) and 
0(3686) produced inside a heavy nucleus will be stud- 
ied via the dilepton spectrum emitted from their decays. 
While the lifetime of J/ has been shown to remain con- 
stant in nuclear matter, those of 0(3686) and 0(3770) 
are reduced to less than 5 fm/c due to increases in their 
width in nuclear matter [33]. Therefore, an appereciable 
fraction of produced excited charmonium states in these 
experiments are expected to decay inside the nucleus [40] , 
leading to an observable shift of the peak positions in 
the dilepton spectrum. The observation of such shifts in 
the masses of excited charmonium states in these experi- 
ments would give us valuable information on the nontriv- 
ial changes of the QCD vacuum in nuclear medium and 
on the origin of masses in QCD. 
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